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Qh Abstract 

"^N The BMSSM framework is an effective theory approach that encapsulates a variety of extensions 

C^ beyond the MSSM with which it shares the same field content. The lightest Higgs mass can be much 

^.^ heavier than in the MSSM without creating a tension with naturalness or requiring superheavy stops. 

i-C The phenomenology of the Higgs sector is at the same time much richer. We critically review the 

I properties of a Higgs with mass around 125 GeV in this model. In particular, we investigate how the 

*r~^ rates in the important inclusive 2j channel, the 27 + 2 jets and the ZZ — >■ 4? (and/or WW) can be 

r^ enhanced or reduced compared to the standard model and what kind of correlations between these rates 

are possible. We consider both a vanilla model where stops have moderate masses and do not mix and 

^N) a model with large mixing and a light stop. We show that in both cases there are scenarios that lead 

to enhancements in these rates at a mass of 125 GeV corresponding either to the lightest Higgs or the 

heaviest CP-even Higgs of the model. In all of these scenarios we study the prospects of finding other 

signatures either of the 125 GeV Higgs or those of the heavier Higgses. In most cases the f r channels 

are the most promising. Exclusion limits from the recent LHC Higgs searches are folded in our analyses 

while the tantalising hints for a Higgs signal at 125 GeV are used as an example of how to constrain the 



o 

^SJ BMSSM and/or direct future searches 
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1 Introduction 

The latest upturn in the hunt for the Standard Model Higgs came at the very end of 2011, when 
the ATLAS and CMS collaborations both presented their analyses with the latest dataset (see 
[TJ[5]) and showed that, besides an exclusion limit that is driving the Higgs mass in a very thin 
region 115-131 GeV, there might be the possibility of a Higgs signal around nih — 125 GeV. If 
these results are confirmed they will mark the crowning of the Standard Model especially that 
this mass range is in excellent agreement with the indirect limit from the global electroweak 
precision measurements that test the inner working of the model, shall we say the theory now, 
at the quantum and renormalisable level. 

With the fact that no new particle outside the SM has been discovered, a lone elementary Higgs 
may bring back the issue of naturalness. A known solution to this problem is supersymmetry. 
No wonder that the first flurry of articles after the announcement of a hint of a Higgs with mass 
around 125 GeV were from aficionados of supersymmetry. 125 GeV is an almost lucky strike for 
its minimal manifestation, the MSSM. Almost lucky, because in the MSSM this value is on the 
heavy side, requiring large values for the mass scale in the stop sector[31IU[51|Sl[71IHl[ni[in]- This 
large scale brings back again the naturalness problem. Extended models of supersymmetry [TTJ[T2] 
fare better from this point of view, in particular one of the simplest versions, namely the next-to- 
minimal version, the NMSSM[T31[Ill[ini[Z], can quite naturally provide a Higgs with mh = 125 
GcV. These more natural extended models allow also a richer phenomenology in the Higgs sector 
than in the MSSM. For example, scrutinising the data suggests that the signal corresponds to 
larger production rates than in the SM in the 27 channel. This is extremely difficult to attain in 
the MSSM [21 [7] where apart from the large values for the scales in the stop sector to obtain the 
Higgs mass, one needs the collaboration of staus with the mass of the lightest stau as low as what 
is permitted by LEP but otherwise very large parameters for other scales in this sector also[5]. 
There is more fiexibility with the NMSSM[13l[7l[T6] when an increase in the diphoton rate comes 
easily as the consequence of a drop in the coupling of the Higgs to bottom quarks due to mixing 
with the singlet component. Although we must stress that, considering the significance of the 
results in the different channels and the two experiments, one should take great care in drawing 
any hasty conclusion, it must nonetheless be admitted that the hint of a signal is at the moment 
still compatible with a SM interpretation. 

One could attempt to better fit the data within the most general effective Lagrangian de- 
scribing the Higgs interactions within an extension of the SM[T71[T51[T^. Fits to the anomalous 
Higgs couplings could then be performed. This could be useful in constraining some underlying 
models but this approach may not be as constraining when the underlying model tightly relates 
the different anomalous couplings of the 125 GeV Higgs that may also have a bearing on other 
(heavier or lighter) Higgses that should be included in the picture. A case in mind is the Higgs 
sector of the BMSSM 23 [SB [H HI [Ml Hi] which is a general effective Lagrangian approach to 
the MSSM first introduced to address the naturalness problem [251 [27] . Masses up to 250 GeV 
can be naturally achieved for the lightest Higgs of the model[2Hl[2Hl[3nilSIl[Slll3S]- Most recent 
analyses of the model taking into account ATLAS and CMS data now only allow masses below 
140 GeV for the lightest Higgs [521155] . Naturally there is no problem for the model to generate 
lightest Higgs masses around 125 GeV. At the same time Higgses with this mass, in this model, 
have a rich phenomenology with properties that can differ from those of the SM and the MSSM. 
It is therefore very important to review the signatures of the different manifestations of these 
Higgses. In particular it is crucial to see which channels, in particular the most sensitive ones: 
77, ZZ{41), 77 H- 2 jets, in this mass range see their rates either enhanced or reduced and study 
correlations between these rates. This is the main purpose of this paper. 



Apart from the profile of the Higgs with ruh — 125 GeV in terms of a signal in different chan- 
nels, we take different constraints in particular non observation of signals of other Higgses with 
the present LHC luminosity. We then take as an example a situation where an enhancement in 
the inclusive diphoton channel and the ZZ — ^ 4/ is confirmed and investigate what consequences 
on other signals either of the same Higgs in other channels or other heavier Higgses are to be 
expected when more data is collected. We also review the role that the stop sector can play 
in the BMSSM. In these models light stops can very easily give mh — 125 GeV, however very 
light stops with large mixings can change in an important way the correlations between the 77, 
ZZ{Al), 77 + 2 jets channels. We also investigate whether the 125 GeV Higgs could correspond 
to the heaviest CP even Higgs. The article is organised as follows. In section [2] we briefly de- 
scribe the BMSSM set up and the main characteristics of the Higgs that ensue, and we define 
the parameter space and how the experimental data is incorporated. In section [3] we consider 
the case where the lightest CP-even Higgs, h, has mass rrih = 125 GeV. In section]!] we turn to 
a scenario where it is the heaviest Higgs, H, that has mass niH = 125 GeV. Section [5] collects 
our conclusions. 

2 BMSSM description 

The BMSSM is an effective theory that builds upon the MSSM by the addition of higher order 
operators. This means that it shares the same field content as the MSSM, in particular in the 
Higgs sector the physical states are the 2 CP-even Higgses (the lightest h and heaviest H), 
the CP-odd A" and the charged Higgs H^ . The higher order operators represent the effect of 
physics beyond the MSSM that has been integrated out and is characterised by a mass scale 
M. The effective theory is then an expansion in powers of 1/M of the Kahler potential K and 
the superpotential W. An exhaustive set of leading order operators has been catalogued in 
|34j but what concerns us here are those operators that have most impact, namely those that 
involve the Higgs sector. Indeed, it is well known that the Higgs sector of the MSSM is very 
much constrained and it is thanks to the radiative corrections that the MSSM has survived so 
far. Likewise perturbing a little through the introduction of these higher operators changes the 
phenomenology of the MSSM quite drastically and improves the naturalness argument. The set 
of leading order operators, beyond the MSSM, is quite restricted [20 l |2 H [22 t [23 l [24 ] 



W,s ^ Cij^ {Hi.H2)\ (1) 

+a^^{H,.H2) (hI.hI) + -^ (a5Hle'''Hi+aeHle''-H2) (h,.H2 + H\.hI) .(2) 



where Hi , H2 are the two Higgs superfields in the gauge basis, assuming the hypercharges to be 
Yi = — 1,^2 = 1- The new parameters Ci: Oi (* = 1--6) are the remnant of the New Physics. 
In order to account for supersymmetry breaking of these extra operators the simplest way is to 
convert those effective coefficients into spurious (see [29]) : 



Ci ^ Cio + Cii^sO^ (3) 

ttt — > fljo + aiimsO"^ + a*imj) + ai2m^J^ d^ . (4) 



We have introduced the spurion breaking mass rtis as a book-keeping quantity so that all effective 
coefficients are dimensionless. One finds that the leading corrections are of order nig/M and 



ji/M , where /i is the usual supersymmetric Higgs mixing parameter. At order 1/Af^ corrections 
are of order {mg/M)'^, (fx/M)'^ and {iinis/M'^). There are also corrections of order w^/Af^, 
V is the Standard Model vacuum expectation value. We assume the underlying theory to be 
approximately supersymmetric so that vris/M is small and the effective expansion well behaved. 
In our analysis we take 






= — = 0.2, M = 1.5TeV. 



M 



(5) 



Note that ^ is (naturally) small, /i 



300 GeV. 



We refer the reader to references [20l EH ESI [33] for the full derivations of the Higgs masses 
and couplings, in particular how the cross sections and branching fractions, either tree-level or 
loop induced, are computed i55]. Here we will just outline the main features. 

With the coefficients (C, a) in the range —1, 1, masses for the lightest Higgs up to 250 GeV can 
be obtained. The largest values are obtained for the smallest tp in the range 2 — 40. tp is the ratio 
of the vacuum expectation values of the Higgs doublets. The coupling VVh/H {V = W^,Z) 
can not be larger than 1 but can be very much reduced. The implementation imposes a custodial 
symmetry, then the WWh and ZZh are to a very good approximation equal, tth is not much 
affected especially compared to bbh that can be either greatly reduced or enhanced. This has 
important consequences for gluon fusion and the branching fraction of the Higgs into photons 
and ZZ in the range of interest around 125 GeV. 

LEP, Tevatron and in particular the LHC data have put dramatic constraints on the model[33] . 
Taking into account the status of the searches at ATLAS and CMS with a dataset up to 2.3 fb"-*^ 
[55j (the parameter space and the experimental constraint being discussed in the following part), 
the lightest higgs, h can not be heavier than 140 GeV, see fig. [l] Note that the data allow for 
the heaviest CP even Higgs, H, to be quite light with mn — 125 GeV, with a light Higgs, h, 
escaping the LEP bounds. There is therefore the possibility that the the hint for an excess at 
125 GeV may be due to H with h being more elusive, a possibility that we will address. 
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Figure 1: We plot here the allowed region in the mH,rnh plane with the 2 fb ^ dataset from 
LHC. For more details, see 



The production channels are also modified. Gluon fusion (gg — > h) can be either greatly 
enhanced or suppressed : when normalized to the SM expectation its value ranges from 0.1 



to higher than 5. Vector boson fusion, VV — >■ h/H, on the other hand can only be reduced. 
The b fusion {bb -^ h) can be greatly enhanced or reduced with respect to the SM prediction. 
It is important to note that the production mechanisms that contribute to the same inclusive 
channel are not rescaled in the same way : hence the differential distributions of the BMSSM 
Higgs production will not be exactly the same as the SM ones since each production mode can 
lead to different kinematic properties. This remark should be kept in mind when exploiting 
SM limits based on inclusive cross section. We discussed this issue at some length and gave 
recommendations in 33 . 

For the decays of the Higges, the branching ratios can also be enhanced or suppressed, 
however these changes are not independent from how the production rates are affected. Since 
we are interested in the mass range around rrih = 125 GeV, we plot in figure [2] the normalised 
value of the gluon fusion as a function of the normalised branching ratio into photons for /i, with 
red (blue) points for R^^ < 1 {R^^ > 1) in the range 122 < rrih < 128 (GeV). In this scenario 
the stop sector has little impact. The ratio R^^, to be better specified later, corresponds to a 
good approximation to the product of the normalised gluon fusion and normalised branching 
ratio into photons. As can be seen, enhancing the branching ratio into photons has a strong 
consequence on gluon fusion. An enhanced branching fraction into photons forces Ogg^h to have 

1.2. 
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an almost SM rate. For the highest value of Br{h — > ^^j)/Br° {h — > 77), agg^hl<^ 

This slight increase , together with the structure of the correlation (the two branches) see fig. l2| 
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Figure 2: We show here the gluon fusion normalised to the SM expectation as a function of 
the normalised branching ratio to photons in the BMSSM, with 122 < to/j < 128 (GeV), before 
applying collider constraints. The effect of the stops is negligible here. 



can be understood as follows. In this set-up gg ^ h is driven mostly by the top contribution 
that is practically standard model like and the bottom contribution that can vary greatly. In 
the SM, the small b contribution interferes destructively with the dominant top contribution. If 



now we parameterise the strength of the hbb in the BMSSM as h^^^^^bb — x^h^^bb, then we 
approximately have 



a 



I _ (1 - efctXfc 



l2 



< ew < 1, (6) 



£64 is the relative contribution of the bottom compared to the top, at the amplitude level. Largest 
values of Br{h -^ ■yy)/Br^^ {h -^ 77) are due to a much reduced h ^^ bb width and hence small 
\xb\. Then for \xb\ < 1 



(^gg^h _ (1 - <^btXb)^ 
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gg^h 
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l + 2e,,t(l-a;b) > 1, 

1 + 2eht Xb -^ 0. (7) 



which is the contribution of the bottom contribution. The bell shape is a reflection of the 
quadratic dependence (1 — Cbtxt)'^- We should keep this characteristic in mind for the rest of the 
analysis. Therefore despite what might seem to be a relative freedom with the new parameters 
introduced by the BMSSM, some predictions are rather constrained. Let us add a comment 
about the reduction of the hbb coupling. With mh ^125 GeV, even a relatively large reduction 
which is accompanied by an almost similar reduction of the hfr coupling, might still not affect 
so drastically the Br{h — > 66) nor the Br{h — )■ tt) since T{h — > 66) might still remain dominant. 

2.1 The parameter space 

tfj and TTij^o (the mass of the pscudoscalar Higgs) are varied in the following range. 

t^e[2,40], m^o e [50,450]. 

Considering the impact of the third family on Higgs physics, we allow some flexibility in the 
stop sector. Two implementations will be considered. 

• Model A: A no mixing scenario, where all soft masses of the third generation squarks are 
set to Mu3j^ = Md3f^ = Mq^ = 400 GeV with no mixing, in particular the stop mixing 
parameter At — 0. For these values the masses mi , mi are around 400 GeV and the 
mixing angle vanishes. This is taken as a standard case, where stops are not too heavy and 
in the set up of the BMSSM their effect is not so important. 

• Model B: A maximal mixing scenario where one of the stop is light ttIj^ ~ 200 GeV. We 
will take to= G [300, 1000] (GeV) and sin20j = —1. The heaviest stop mass is taken as a 
free parameter. This will have important consequences in the production of the Higgses 
and their decays. Note that, in a generic model, a 200 GeV stop can still escape all current 
collider limits. 

For the rest of the MSSM parameter space, the following values have been used : all soft 
masses are set to Mgait = 1 TeV (except for the third generation), ^ and M^ (the SU{2) gaugino 
mass) are set to 300 GeV, Mi (the U{\) bino mass) is fixed by the universal gaugino mass relation 
Ml = I tan2 ewM2 ~ M2/2, and M3 = 800 (the 5C/(3) gaugino mass) GeV, cos^ Ow = M'^/Ml. 
All trilinear couplings are set to (except At). 

For the purpose of this paper this last set of parameters has practically no impact, so we 
could have taken lower values of Mgoft and M3 compared to the new scale M — 1.5 TeV. In 



scenario B, the largest value of the heaviest stop m^^ = 1 TeV that we allow in the scan should 
be regarded an extreme example, not only from the point of view of naturalness but also because 
it is not far heavier than the new scale M = 1.5 TeV. This said we must emphasise that most of 
our study is done with nii — 600 GeV. Furthermore the heavy scale M can be enhanced with 
little change to our results provided one keeps fixed the ratios ras/M^ji/M . 

The effective coefficients, (C, a) will be varied in the range [—1, 1]. The following constraints 
were applied on each point, in order to constrain the effective parameter space (see ^3 for more 
details) : 

• Perturbativity check. We check that 1/M^ contributions to rah are small enough. Indeed 
one must make sure that the perturbative expansion in 1/Af is under control. This means 
that on top of all experimental constraints that we take into account, we also check that our 
points exhibit a correct effective expansion by keeping only points where the third order 
{i.e. 1/M^ terms) contribution to the light Higgs mass is within 10% of its value. 

• Electroweak Precision tests. Since the effective coefficients have a non vanishing contri- 
bution to oblique parameters, we verify the consistency with the electroweak precision 
measurement s . 

• LEP and Tevatron Higgs searches. This applies to all Higgses : h,H,A'^,H'^, including 
top decays to H^ . 

• ATLAS and CMS Higgs searches. The list of channels that we take into account and how 
these are exploited within the BMSSM is detailed in[53]. As concerns the exclusion limits 
exploiting the tt channels we had also included in[33j the MSSM analysis of ATLAS and 
CMS. Concerning the mass range 122 — 128 GeV and an interpretation in terms of a signal, 
we detail our approach in section |2.2[ 



In this work we do not include constraints based on flavour physics {{g — 2)^, Bg — )■ fj,'^ fj,~ ,..) 
and dark matter (relic density) constraints. For the former, this will introduce some extra model 
dependence from extra operators in the BMSSM, outside the Higgs sector. For the relic density it 
is known that the prediction can change drastically if we change the cosmological model. These 
issues may obscure the conclusions on the interpretation of the possible Higgs signal. We will 
certainly come back to these effects in future investigations. 

2.2 Input from the LHC 

In order to use the results from the ATLAS and CMS collaborations, we have used the following 
ratios 

p g"pp-»/i-»XX 1 ncxclusion <^pp^H^XX /gx 

^XX - -SM and KxX = ^excluded 95% ' ^^1 

where XX denotes a particular final state (say the inclusive 27) . (jcxciuded 95% g^^nds for the 95% 
C.L. excluded cross-section reported by the collaborations. In practice the Rxx will be used 
in the signal case, to compare with the best fit /t - of the so called signal strength fi- given by 
the experiments. In eq. [8j ft. in the BMSSM will refer either to the lightest or heaviest CP-even 
Higgs. i?3fx"*''°" ^ill b^ used in the no-signal case as a measure of the sensitivity of the search, 
here H stands for all Higgses not contributing to a signal in the mass range 122 — 128 GeV. It 
will also be shown as a measure of the luminosity needed to see the effect of a particular Higgs 
in a certain channel in the future. For Rxx the most important channels are the inclusive 27, 
ZZ -^ 41 (ATLAS and CMS) and the exclusive 27 -I- 2 jets (CMS). For R''^f'^'°'^ we use the full 



dataset across all Higgs masses covered by the experiments with the current collected luminosity 
of 4.9fb~^. As we said i?3fx"^'°" gives a measure of the luminosity needed to uncover a signal in 
a new channel. The current i?3fx"*''°" i^ based on a collected luminosity of about 4.9fb^ for the 
LHC running at 7 TeV. To extract from the plots for /JJ^^"*""" that we will show the approximate 
luminosity that will be needed to uncover a potential signal, one can base the rough estimate 
on a simple rescaling of the luminosity. For example in the plots we will show, a channel with 
i?™x"'''°" = 0.1 will require a luminosity of about 500fb~^ to be observed. With a luminosity of 
30fb~ only those channels with i?3fx"'''°" > 0-4 may be accessible. 

We decided not to carry out a fit of the model to the data (as was done in [TOl [TH] among 
others) for the reason that so far the signals from both collaborations are not so easy to reconcile 
and also because in view of the quality of the results this exercise is far too premature. We chose 
instead to focus on assessing the limits on the flexibility of the BMSSM with a Higgs in the 
range 122 < rrih < 128 GeV, a range broad enough that it takes into account the hints from both 
experiments including the uncertainty on the mass in each experiment. We leave a complete 
computation of the compatibility of the model with the data for the future, hopefully more 
precise, set of data. By using those ratios instead of computing the complete likelihood function 
(as was tried in |18|). we have implicitly assumed some approximations that are discussed in |33j . 
In the inclusive channels, we neglect for instance the effect of acceptance and efficiency cuts due 
to the change in the relative contribution of the different production modes (gluon fusion, vector 
boson fusion and associated production, b quark fusion) compared to the Standard Model. In 
other words, we are taking ratios of inclusive cross-sections. For analyses that are truly inclusive, 
such as the h — >■ ZZ — ^ 41 there is hardly any difference (see |53]), but the situation changes when 
channels are divided in subchannels with different final states. The h/H — )■ WW + 0/1/2 jets 
comes to mind. The difficulty is that each subchannel that contributes to the inclusive cross 
section has a different efficiency and the fact that most often the rates for the different channels 
in the BMSSM are not rescaled by the same factor. The reason for these approximations is that 
most of the experimental quantities, such as the efficiency and the acceptance for each production 
mode or the cross-section limits for each subchannels, are so far unavailable, forbidding hence an 
exact interpretation. So we reiterate our recommendation that efficiencies be provided (see the 
Les Houches Recommendations for the Presentation of LHC results [S^ for a detailed discussion 
on the subject.). However, for the time being those approximations do not prevent a qualitative 
survey of the would be Higgs signal properties. 

Having said that, CMSjSTj has most recently provided data for a more exclusive observable, 
27 -|- 2jets. The latter is more sensitive to the production through vector boson fusion, with 
the Higgs subsequently decaying into 2 photons. We take this observable into account. The 
contribution of the gg —> h to this channel is however not negligible. Fortunately the CMS 
collaboration fST does provide the overall acceptance/efficiency (product of the two) for each 
of the important channels (this applies also for the BMSSM), albeit for a single Higgs mass, 
TO/i = 120 GeV. The overall acceptance is quoted as 0.15 for production through Vector Boson 
Fusion (VBF) and 0.005 for gluon fusion. We therefore simulate the ratio Ry^+2 jots as 

_ 0.15 (TvBF + 0.005 (Tgg^h BR^^ 

"^+'^'*^" 0.15 a^^^p -I- 0.005 a^^/f,, X ^^^ ^yj 

We checked that this parametrisation of (7-y^+2 jots when folded in with the SM cross sections for 
the LHC at 7 TeV [55] and taking into account the luminosity quoted by CMS reproduced quite 
exactly the number of selected events given by CMS [37] . We assume that this parametrisation 
that was verified to be excellent for nih = 120 GeV still holds to a very good degree in the range 
122 <mh< 128 GeV. 

In order to determine if a Higgs is excluded, we use a quadrature combination of all /^exclusion 



from all channels of both experiments and checked if the result was below 1. We apply the 
test on all Higgses {h,H,A'^,H^), adding ratios when the Higgs bosons are degenerate, before 
determining whether the parameter point is allowed. 

The signal condition we will require is the following : first we require a Higgs (heavy or light) 
boson in the range 122 — 128 GeV. Then, we discard any exclusion limit on this Higgs in the 77, 
ZZ and WW channels (since we want it to be a signal in those channels). We apply on all other 
channels and on all other Higgses the latest exclusion limits, obtained with the 4.9 fb~^ dataset. 



2.3 Signal features, data 

The data that is most indicative of a possible signal is the following (uncertainties correspond to 
the Icr band) 

• ATLAS [39]: 

The ATLAS collaboration records a combined (all channels) signal strength of 1.5 l^'g at 
iTih = 126 GeV. It may be considered as most revealing in channels with best resolutions 
on the Higgs mass: 

• The inclusive 77 channel where the signal strength is 2_Qg 

• ZZ — )■ 41 channel where the signal strengths is 1.2^Qg compatible with WW — > llvi/, 
though the WW channel has a worse mass resolution. 

• CMS collaboration reports a combined signal strength of 1.2^q'4 at nih = 124 GeV 



• 



In the 77 channels, the first CMS release with 4.9 fb^^ was based on an analysis with 
four subchannels that gave a signal strength of 1.7 ± 0.8 at m/j — 123.5 GeV (see ref 
([lO]))- The updated release added a dijet-tagged subchannel 77+2 jets yielding by 



• 



itself a signal strength of 3.8_]^g. The combination of the five subchannels yield a 
signal strength of 2.11^7 (EZ]). 

For the ZZ — > 41, 0.5^q ". Note that the mean is low, moreover the mean value for 
771^ is at 126 GeV. 



• the bb and t^t channels analysed by CMS[40] in the mass range 122 — 128 GeV have 
so much uncertainty that they are of little use in the present analysis. 

Let us take breath and emphasise again that there is still much uncertainty in these results, 
some of which may not help in drawing a coherent picture, expect perhaps in the 77 channel. The 
signal strengths are compatible with a Standard Model Higgs, however it is tempting and in any 
case educative to entertain the idea that some non standard Higgs scenario is emerging. What 
is very interesting is that the different channels and subchannels will allow, when measured with 
better precision, to discriminate between different models and implementations of the BMSSM. 
Most probably a first step in this discrimination in this mass range will be performed with 
77, VV, 77 + 2jets perhaps also with the incorporation of the tt channel. In the case of a multi- 
Higgs system this will be done in parallel with searches for other Higgses. In the rest of the 
paper we will investigate what kind of correlations between these observables are possible within 
the BMSSM, for example whether enhancements in all channels are possible. 




Figure 3: Allowed region in the plane ruh, Rj-y. The blue line represents the ATLAS best fit for 
the signal strength, and the dotted lines are the one sigma deviations from this value in model A. 



3 h as a. signal in the 122 — 128 GeV window 

3.1 Model A: No light stops, no mixing in the stop sector 

Fig. |3] shows that with the current data, the BMSSM yields a production rate in the inclusive 
pp — > /i — > 77 that can be quite small (as small as 0.1), and hence unobservable with the current 
luminosity or in the very near future. More interestingly there is however no difficulty finding a 
signal in this channel that is up to 3.5 times that of the SM. There is a very strong correlation 
with the signatures in the other promising channels: VV = ZZ — >■ Al and the 27 + 2 j'eis, see 
fig 111 With small differences we have -R^-y ~ Rzz ^ ^77+2 jots- Rates above those of the SM are 
mostly driven by reduction in the width of to bb which increases all channels. This is trivially seen 
for the 27 versus ZZ channel. In the case of the j"f/"fj + 2 jets correlation, when the rates are 
above those of the SM, the inclusive channel is higher by 20% or so, this (and the appearance of 
two branches) rests on the same argument that we put forth in section p] about the contribution 

= 3 is very 



>4; 



X , JL.-^^ — ^ , jX- 



77 + 2 jets 



of the b quarks. Therefore a configuration with Rzz- 
much disfavoured in Model A. 

It is important to stress that the characteristics we find in these scenarios occur for all values 
tfi , even if statistically, with a simple scan, the population with smaller tp is larger. 



3.2 Model B: a light stop with large mixing in the stop sector 

It has been known for some time[lTl |31] that, within the MSSM, light stops endowed with a 
large mixing can drastically reduce the gg induced production. Even if this is accompanied by an 
increase in the decay width to photons, the combined effect can be a large drop in gg — ;> ft, — >■ 77. 
This effect is encapsulated in the coupling of the stops to the Higgs. The coupling of the lightest 
stop, ti, gfiiiii writes in the large ruj^o limit 
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(10) 



9i is the mixing angle of the stops. The ^2 coupling is obtained through ti -i-> ^2- The non mixing 
term m^ adds up with the top contribution, whereas the mixing term interferes destructively 
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Figure 4: Le/t panel: correlation between R^^, Rzz and R-y-y+2 jets for 122 < mh < 128 GeV. 
Right panel: Imposing R^^ h = 2.0 ± 10% (points in red) and R^^ h = 2.0 ± 1% (points in black) 
we show the correlation in the plane Rzz o,nd R^^^2 jets- Both figures are for model A. 

with the top. For large mixing with large enough gap between the two stops masses this means 
that a reduction in gg ^ h occurs but accompanied with a more modest increase in the h — > 77 
due to the W loop. Of course the Br{h — >■ 77) can be much more efhciently increased if a drop 
in h —i' bb occurs as within the BMSSM. Therefore we see that by letting light stops jump into 
the game and keeping a ratio in the 77 channel higher than the standard model, the correlations 
between the different channels will change. 
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Figure 5: Allowed region in the plane ruh, R-y-y. The blue line represents the ATLAS best fit for 
the signal strength, and the dotted lines are the one sigma deviations from this value in model B 
with maximal mixing and with mi^ = 600 GeV. 



m^^ = 600 



We first note, see fig. pi that in the maximal mixing case siv?{26i) — 1 and with 
GeV, R^^ is reduced sornenow compared to the non mixing case of model A, however one still 
obtains enhancements of a factor 2 (and more) compared to the SM. However, now the 77 + 2 jets 
can be much higher than the 77 channel, whereas previously we had i?^-y+2 jets 
now for the same value of R^^ R^'f+2 jots = 2.5, see fig.l6j Moreover the weight between R 



1.5 for Rj^ — 2, 



and Rzz has been inverted, we now have R. 



■77+2 jets 



77+2 jets 

> Rzz- Scanning over mf from 300 GeV 
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to 1 TeV will open up more possibilities for the correlations between these channels. The results 
of this scan are shown in the right panel of fig. [6] For example imposing that R^^ = 2.0 ± 10% 
one can obtain R^-^ + 2 jots = 3.8 together with Rzz — 1-3. We can therefore recover values that 
correspond to the best fits for these observables obtained by the two collaborations. We stress 
again that this is illustrative and shows how much flexibility in the model can be introduced. 
While in the case of no-mixing in the stop sector all channels seemed to have nearly trivial 
correlations, raising the mixing will in most cases raise the 77 + 2 jets channel compared to the 
77 channel, and also decrease the ZZ — )• Al channel with respect to the 77 one. 
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Figure 6: Left panel: correlations between R^-^, Rzz o-nd R'yy+2 jets for 122 < rrih < 128 GeV 
in the maximal mixing scenario of model B with ruf^ ~ 600 Ge V. Right panel is a subset after 
imposing R^^ ^ — 2.0 ± 10% (points in red) and R^^ h = 2.0 ± 1% (points in black) in the plane 
Rzz cind R^-y+2 jets iiT' rnodel B scanning in the range m^ G [300, 1000] (GeV) 



3.3 Prospects for other signals 

Although an unambiguous signal refuting the SM would be, in the case where the signal at 
mil — 125 GeV were confirmed, a precise determination of the signal strength above (or below) 
the SM expectation, such a precision may require some time. At the same time as the luminosity 
increases other channels and signatures may become sensitive in corroborating the signals with 
ruh ~ 125 GeV. These channels could either be other channels where the same Higgs with mass 
125 GeV takes part or channels affecting the other Higgses of the model. In the first case, the 
other allowed decay modes are tt and bb final state, however if the trend towards an increase in 
the 27, ZZ and 27 + 2jets is reinforced requiring a reduced hbb (and consequently Htt) in the 
BMSSM, the tt and bb whose current sensitivity in the SM is quite low will require substantial 
increase in the luminosity. 

To pursue this investigation about the prospects of signals in other channels, we keep for the 
sake of illustration those models compatible with 



1.2 < R^^ < 2.9 



k 0.5 < Rzz < 2.4, (11) 

which is the one sigma band obtained by the ATLAS collaboration and show the different 



pexclusion 



Again, eq. 
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is an arbitrary choice, taken for the sake of concreteness. One should 
keep in mind that as more data is collected, this requirement will become stronger or perhaps 
even totally irrelevant. 
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3.3.1 Non-mixing scenario 
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Figure 7: Discovery perspective for the channel h ^ tt . 
SM. 



The line is black corresponds to the 



We start with the tt signal of h. We see in fig. ul that /^c'^ciusion jg always below 0.4, with 
a concentration below 0.3 that corresponds also to the expectation from a SM Higgs, therefore 
a luminosity in excess of 30fb~^ is needed in the most favourable cases. Most cases will require 
much more luminosity, up to 500 fb""'^ in the worst case. Incidentally we note that this channel, 
despite the reduced hrr coupling, can be above that of the SM, which shows that a reduced Htt 
does not mean a large drop in the tt branching ratio, moreover the production cross section can 
be larger than in the SM. 
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Figure 8: Discovery perspective in the tt channel through the heavier Higgses: in the case where 
A^ and H are degenerate within 10 GeV (left pannel) and when A^ and H are not degenerate 
(right panel). In the latter the correlations between the two signal is shown. In the panel on the 
left, the SM case is shown in black. The different shades for the BMSSM correspond (from left 
to right) to cases with tp < 10 (red) , 10 < t^ < 20 (green) and tp > 20 (blue). 

Would the other Higgses be more sensitive? The answer can be drawn from fig. [8j Some sce- 
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narios can be probed with little increase in the present luminosity. Generically, high tp [tfj > 20) 
will be probed within the next 30 fb^-'^, while low tp {tp < 10) could be quite hopeless if the 
heavier Higgses are heavier than 400 GeV. We find that /j<=xciusion ^ q g ^^^ reached in cases 
where A° and H are close enough in mass to be degenerate {\mjio — mnl < 10 GeV), yielding 
thus a single signal, /^exclusion ^ q g jg j-g^ched also when the degeneracy is lifted, in which case 
one expects both signals to be revealed with roughly the same luminosity, see the correlation in 
fig.^ Models with niH ^ m^o < 250 GeV (degenerate case) show a ratio i^exciusion ^ q^^^ which 
means that the region where the decoupling is not complete between light and heavy Higgses 
could be probed with about 30 fb^^ . In the non-degenerate case, there is of course a loss of 
a factor two, but there is still a lower limit to the exclusion ratio in this mass range. But in 
many models, /^exclusion ^ q^^ ^o g^j^^j jj .^jjj g^ undetected even with a luminosity in excess 
of 30fb~^. This discussion shows that studying the r channel in Higgs physics is crucial. Not 
only it can deliver new signals but can give important information on the parameters of the model. 
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Figure 9: Discovery perspective for other signals : H — > ZZ on the left and H -^ WW on the 
right. No mixing or light stops are assumed. The curve in black is the SM Higgs hypothesis. For 
WW , the curve is out of the hounding box, this confirms that for Higgs masses above the WW 
threshold this channel is very constraining. 

Other channels offer little prospects, apart if Mh ~ 180 GeV where the search sensitivity in 
the clean WW and somehow also the ZZ channel is high, despite the fact that the HWW is 
quite small, see fig. [9J 



3.3.2 Maximal mixing 



In the maximal mixing case, with m^ — 600 GeV, there are few differences. The drop in gg — > h 
is the reason behind the drop in s ensitivity. Subsequently the tt channel of h will be even less 
sensitive, as can be seen in fig. 
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^exclusion ig j^Q^ ^3glQ^ Q 2, 



As concerns the heavier Higgses, the changes are marginal compared to the no-mixing case. 
The best prospects are in the r channels and in the WW channel if mn ~ 180 GeV. The 
corresponding figures are similar to those shown for the no mixing case and we therefore do not 
display them here. 
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Figure 10: Discovery perspective for the channel h 
line in black represents the SM. 



TT in the maximal mixing scenario. The 



4 iif as a signal in the 122 — 128 GeV window 

As fig. [l] makes clear, the BMSSM is compatible with a scenario where it is the heavier of the 
two CP even Higgses, H , which is in the range 122 — 128 GeV and may thus be responsible for 
a signal, while the lightest Higgs h has so far gone undetected. Such possibility, even though 
restrained, has also been evoked in the case of the MSSM[5]. We review such a possibility in 
the case of the BMSSM both in a scenario with no stop mixing and a scenario with large stop 
mixing and light stops. 

4.1 Model A: No light stops, no mixing in the stop sector 
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Figure 11: We show here the allowed region in the plane mu, R^-y (left panel) and the associated 
correlations between R^^, Rzz arid R^-y+2 jets for 122 < wlh < 128 GeV (right panel) in the 
scenario with no-mixing . 

The statement we have just made can be made more quantitative. Solutions with 122 < 
niff < 128 GeV correspond to a situation where all three Higgses are light in the sense of being 
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all three below the WW threshold, m^ < 120 GeV m^o < 160 GeV. We find that some features, 

for the signal observables, are to a large extent similar to what we have found in the case of h. 

In a way the h and H have swapped their role as to which is the SM-like, SM-like as concerns 

the VVH/h strength. Indeed, this is illustrated in fig. 

as 3.5, there are correlations between R^-y, Rzz and R 

cases, but not all as was the case for 122 < m/j < 128 GcV. In this case, there is some spread in 

the correlations between Rzz and R-y-y+2 jots, see fig. 
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R^ 



-y^ can still reach values as large 
77+2 jets with Rzz > R'fy+2 jots in most 
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Let us turn to other characteristics of these BMSSM scenarios and how they could show 
up in other observables. Another mode where H could be observed is the H ^>- tt channel. 
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Figure 12: jiExciusion j^^ ^22 < iriH < 128 GeV. No mixing is the stop sector is assumed. The 
line in black is the SM. 

Fig. [12] suggests that prospects here might be better than for h giving a signal in the range 
122 <mh < 128 GeV. Indeed, there are solutions with RE^dusion ^ q 5 ^j^^^^ would need about 
20fb~^ to be uncovered. 

Observability of the other Higgses shows, in many cases, very good prospects, gain in the tt 
channels, j^Exciusion ^ q g ^^^^ obtained, see fig. [Ts^ . Therefore it is worth pursuing searches of 
/i, for ruh < 120 GeV in the tt channel. A^ could also be uncovered with the same luminosity, in 
fact fig. 13 shows the correlation between h and A'^ in the tt channel. There, of course, remains 



also many situations with R 



Exclusion 



< 0.2 that would be difficult to decipher. 



4.2 Model B: a light stop with large mixing in the stop sector 

We now turn to the maximal mixing case and restrict ourselves to rrij:^ — 600 GeV {m^ — 200 
GeV). Compared to the previous case without mixing one notes that there is a reduction in 
RytH- This is mainly driven by the drop in gg — > H , see fig. 14 very low values are also due to 
quite small Htt couplings. 



The most noticeable change is the correlation between Rzz and i?^^+2jets, see fig. 15 We 
now easily find i?^-Y+2 jets > Rzz- The spread in this correlation has increased. One can find 



scenarios with Rzz < 1 even for R^^ > 2. For R^^ ^ 2, R, 



^77+2 jets 



> 2 is attained. 



The visibility of A'^ and h is little affected by the stop mixing. Our conclusions are little 
changed. Again it is very important to pursue the search in the tt channel. 
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Figure 13: Ji±xciusion ^^ ^ function of mh for 122 < iriH < 128 GeV in the case where A" and 
hare not degenerate within 10 GeV (left panel). The right panel shows the correlations between 
the h and A'^ in the t channels. No mixing in the stop sector is applied. The line in black in the 
left panel represents the SM. 
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Figure 14: RttH iri the same mass 122 < m^j < 128 GeV in the case of maximal stop mixing 



and ' 



600 GeV. The line in black is the SM. 



5 Conclusion 

Despite the fact that no sign of supersymmetry has been found so far, the BMSSM framework is 
a very efficient set up that extends the realm of the MSSM in a most natural way as concerns the 
reahsation of the Higgs. In the MSSM framework there is some tension with naturalness for a 
Higgs mass of 125 GeV that requires heavy stops, in the BMSSM this is not an issue. Although 
one must exercise extreme care with the so called tantalising hints for a Higgs signal around this 
mass, 125GeV, it is extremely important to scrutinise the properties of the Higgs with such a 
mass in many models, in particular the BMSSM which represents an effective implementation of 
a variety of supersymmetric models having the same field content as the MSSM. The tantalising 
hints have come also with the temptation, even if premature, of reading from the results of 
ATLAS and CMS, despite the uncertainty of the measurements, that the signals in the inclusive 
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Figure 15: We show here the allowed region in the plane mn, R-yy (left panel) and the associated 



correlations between R^^ 
maximal mixing and mf 



Rzz and iJ^^+a jets 
= 600 GeV. 



for 122 < niH < 128 GeV (right panel) with 



27 channel, the 27 + jets and perhaps in the ZZ — > Al to be higher than what is expected from 
the SM. Such scenarios are practically impossible to attain in the MSSM, a possibility that has 
been entertained would make the naturalness argument even more excruciating. It is therefore 
very important to find out whether some configurations, especially those leading to enhancements 
in these most important channels can be realised in the BMSSM. As important is to find out 
how these enhancements or signals are correlated and how different kinds of correlations can be 
realised. We have shown that a vanilla BMSSM where stops are at very moderate masses and 
with little mixing easily allows enhancements in all these channels for mh ~ 125 GeV with the 
constraint that the rate ZZ — > 41 would generally be higher than the rate 77 +2jets. A light stop 
with large mixings in the stop sector offers more possibilities especially as concerns correlations 
between these three important channels. Our study also reveals that although it is easier to have 
such realisations work for the lightest Higgs of the BMSSM, solutions where it is the heaviest 
Higgs that has a mass around 125 GeV also exist. Once a signal at 125 GeV has been confirmed 
a better measurement of the rates, in particular the 27 (inclusive and exclusive) as well as the 
4:1 would narrow considerably the parameter space of the BMSSM. At the same time as more 
precision is achieved and more luminosity is gathered one can constrain the models through 
the other Higgses (those outside the 125 GeV window) but also through other channels of the 
Higgs at 125 GeV. Our study reveals that in both cases it is crucially important and telling to 
investigate the tt channel. We have not folded in the possible constraints from flavour physics 
and dark matter as we have argued that this introduces some model dependencies (including 
those from cosmology) but it is clear now that we have entered a fascinating era. The study 
we have conducted is an example which shows that even before any new direct signal of New 
Physics is discovered, the study of the Higgs, once confirmed, will give important clues on the 
New Physics. We eagerly await more data and analyses from the experiments and we urge, once 
more, our colleagues to provide as much information as possible on the data so that one can gain 
access to the different individual subchannels that make up an inclusive channel. 
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A Addendum: Tevatron and the bb channel 

While this work was being finahsed, the Tevatron Collaborations released new analyses [13] point- 
ing out to a possible signal in VH — >■ Vbb channel with a rate that could be compatible with the 
Standard Model expectation and with a mass that could correspond to where the excesses are 
seen at the LHC. This would seem at first sight to disfavour a scenario where g^f,;, is very much 
reduced. However, one must keep in mind that since the decay H ^ bb dominates for mh — 125 
GeV, a suppression of the coupling by a factor two does not imply a suppression of the branching 
ratio by a factor two. The suppression in much more modest and there can still be a significant 
enhancement of the diphoton channel without suppressing too much the VH — >■ Vbb channel. It 
must be stressed that a more precise measurement of the latter process would really be helpful. 
Indeed, there exists also a correlation between the diphoton (inclusive) channel and this channel, 
as shown in fig [16] 




Figure 16: Correlation between diphoton channel (R-y-y) and the VH — > Vbb (Rii,) in the nan 
mixing scenario. 
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